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An approximate functional equation for effective conductivity a^f f of systems 
with a finite maximal scale of inhomogeneities is deduced. An exact solution 
of this equation is found and its physical meaning is discussed. A two-phase 
randomly inhomogeneous model is constructed by a hierarchical method and its 
effective conductivity at arbitrary phase concentrations is found in the mean- 
field-like approximation. These expressions satisfy all necessary symmetries, 
reproduce the known formulas for a^f f in weakly inhomogeneous case and 
coincide with two recently found partial solutions of the duality relation. It 
means that a^f f even of the two-phase randomly inhomogeneous system may 
be a nonuniversal function and can depend on some details of the structure of 
the inhomogeneous regions. The percolation limit is briefly discussed. 



PACS: 73.61.-r, 75.70.Ak 

The electrical transport properties of the classical inhomogeneous systems 
have an important practical interest. The simplest problem in this region is a 
finding of the effective conductivity erg/ / of an isotropic inhomogeneous (ran- 
domly or regularly) two-phase system, which is a mixture of two phases with 
different conductivities 17^(1 = 1,2). Despite of its relative simplicity only a 
few general results have been obtained so far. In case of weakly inhomogeneous 
isotropic medium pP 

-eff = (a) 1^1 - j . (1) 

For two-phase system (cr) = xai + (1 — x)a2, (cr^) — (cr)^ = 4a;(l — x){a-)^, 
where a; is a concentration of the first phase, (t_ — {ai — fT2)/2 , and (1) takes 
the form 

aeff - a+ (1 + 2ez - (1 - 4e2)zV2) , (1') 

where a+ = (ci -I- cr2)/2, z — cr-/a+ (—1 < 2: < 1), e = a; — 1/2. 

Another general formula is the dilute limit of the Maxwell-Garnett formula 

m 

aeff=ai{l-2{l~x)z), (2) 

where 1 — x <^ 1 is a small concentration of the second phase and a round form 
of the inclusions of this phase is suggested. 
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Then, the exact Keller - Dykhne formula for (Te// in systems with equal 
concentrations of both phases and the exact dual relation, connecting effective 
conductivities at adjoint concentrations x and 1 — x, have been obtained Pl E] 

f7e//(a;, Cri,(T2)o-e//(l - X,ai,a2) = Cri(T2 (3) 

The exact Keller - Dykhne formula follows from (3) at x = Xc = 1/2 

(Teff{xc,(Ji,(72) = ^/o\a^■ (4) 

This formula is very simple and universal, because it does not depend on in- 
homogeneous structure and takes place for regular inhomogeneous two-phase 
systems 0] |S] as well as for slightly nonregular inhomogeneous systems ■ 

Of course, a formula for the effective conductivity at arbitrary phase concen- 
trations has the main interest in this problem. One such approximate formula 
for cTe// has been obtained many years ago in the so called effective medium 
(EM) approximation |H] 

(Te//(e, {a}) = a+ (2ez +^l-z^+ 4(ez)2) . (5) 

where {a} — (cri,cr2)- Though it corresponds to the weakly inhomogeneous 
case it turns out to be a good approximation, when ^ 

It was shown recently that the duality relation (3) together with boundary 
conditions and some assumptions about the properties of Cej^j allow to find 
an explicit form of (Jeff at arbitrary x and two such expressions were found 
P", "5^. In this paper we will represent two randomly inhomogeneous models, 
having their effective conductivities just of these two forms. Another impor- 
tant question appears naturally in this problem: is a formula for the effective 
conductivity of randomly inhomogeneous two-phase systems universal or it can 
depend on the structure of the inhomogeneities of the system? The strong ar- 
guments in favor of the nonuniversality of the effective conductivity in case of 
three-phase regular inhomogeneous systems were derived in the paper |10j . Our 
results for CTg// demonstrate that a formula for the effective conductivity may 
be nonuniversal even in the two-phase case. At the end of the paper we will 
briefly discuss some peculiarities of the percolation limit. 

We start our investigation with a general discussion of the averaging proce- 
dure for obtaining aeff ■ It is easy to see that the effective conductivity will 
depend on a scale I of a region over which an averaging is done. This takes place 
due to the possible existence of different characteristic scales in the inhomoge- 
neous medium. In the most general case there will be a whole spectrum of these 
characteristic scales. This spectrum can be very different: from discrete finite 
till continuous infinite, and is defined by the structure of the inhomogeneities 
of the system. For this reason it can depend on the phase concentration x. 
Suppose, for the simplicity, that the randomly inhomogeneous structure of our 
system has the scale spectrum with a maximal scale lm{x), which is finite for 
all X in the region 1 > x > 1/2 (or 1 — x in the region 0<1 — a;<l/2). Let 
us assume that we know an exact formula for a^f / {x, {a}) of this system, which 
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(b) 

Fig.l. a) An elementary square of the model with round inclusions, b) a lattice of the 
model; the numbers 1,2 denotes squares with the corresponding concentrations. 



is applicable from scales I > Im- It means that this formula for aeff{x,{a}) 
takes place after the averaging over regions with a mean size I > Im and does 
not change for all larger scales I ^ Im- Now consider a square lattice with the 
squares of length II ^ Im and suppose that they have the effective conduc- 
tivities corresponding to different values of the concentrations Xi and X2 with 
equal probabilities p = 1/2 (see Fig.l) 

After the averaging over the scales I > II one must obtain on much larger 
scales I ^ II the same effective conductivity, but corresponding to another 
concentration x = {xi + X2)/2. This is possible due to the similar random 
structure of different squares and due to the conjectures that (1) in this model 
there are only two important scales: a maximum of the maximal characteristic 
scales (we suppose here that Zi ~ Z2 ~ lm{x) , where lm{xi) = k {i = 1,2) ) 
and the lattice square size II ^ h , (2) the averaging procedures over these 
scales do not correlate (or weakly correlate) between themselves. Thus, for a 
compatibility, all concentrations must be out of small region around critical 
concentration Xc, whore li or /,„ (./■) can be very large. Wc will call further 
this set of the conjectures the finite maximal scale averaging approximation 
(FMSA approximation). It can be considered as some nontrivial modification 
of the EM approximation and can be implemented for systems with compact 
inhomogeneous inclusions with finite Im ■ 

From the other side the effective conductivity on scales I ^ II must be 
determined by the universal Keller - Dykhne formula (4). Thus we obtain the 
next functional equation for the effective conductivity, connecting aeff{x, {a}) 
at different concentrations. 



(Jefs{x,{a}) = ^Jae^^ix^,{a})aeff{x2,{a}), x={xi+X2)/2. (6) 
It must be supplemented by the boundary conditions 

0"e//(l, {O-}) = C71, C7e//(0, {a}) = f72, (6') 

The equation (6) can be considered as a generalization of the duality relation 
(3), the latter being a particular case of (6) at + a;2 = 1. It follows from 
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(6) that for z ^ \ and due to the exactness of the duality relation it really 
works at all concentrations x except maybe of small region near x > Xc and 
(72 = (z = 1) (see below a discussion of the percolation limit). It is easy to see 
also that the approximate formula (2) for cre//(a;, {cr}) satisfies equation (6). 
Moreover, one can find an exact solution of this equation. It has an exponential 
form with a linear function of x (or e ) 

o-e//(x, {cr}) = CTi exp(aa; + 6), (7) 

where the constants a, h can be determined from the boundary conditions 

a — ~b, exp 6 = (72 /(Ti . (7') 

Substituting these coefficients into (7) one obtains 

aeffix,{a}) = a, [a^la^f-''^ = o^^a^^-^ (8) 

The solution (8) satisfies all required symmetry relations and exactly coincides 
with the case (a) from [S) |3| ■ 

It is interesting to note that the form (8) means that in the FMSA approx- 
imation one has effectively the averaging of In cr since it can be represented 
as 

IntTe// = (Inf) = a;ln(Ti + (1 — x) ln(T2. (9) 

This was noted already in 3 for the case of equal concentrations x = 1/2. The 
analogous result has been obtained later for random system in the theory of 
two-dimensional weak localization 11 . One can check that (8) reproduces in 
the weakly inhomogeneous limit the universal Landau - Lifshitz expression (1). 
In the low concentration limit of the second phase it gives 

1 — z 

cre//(a;,z) =cri(l + (l-a;) log— — + ...), 1 - x < 1, (10) 

J. ~|~ 2, 

what coincides with (2) in the weakly inhomogeneous case. Note that the ex- 
pansion (10) contains the coefficients diverging in the limit \z\ —f 1. Such be- 
haviour of the coefficients denote the existence of a possible singularity in this 
limit, where the FMSA approximation is not applicable (see below a discussion 
of this percolation limit). 

Now we will construct a hierarchical model of flat isotropic randomly inho- 
mogeneous two-phase system, using the composite method introduced above, 
and find its effective conductivity aeff{x,{(T}). 

Let us consider a simple square lattice with the squares consisting of a 
random layered mixture of two conducting phase with constant conductivities 
CTi {i = 1,2) and the corresponding concentrations x and 1 — x. A schematic 
picture of such square is given in Fig. 2. 

The layered structure of the squares means that the squares have some pre- 
ferred direction, for example along the layers. Let us suppose that the directions 
of different squares are randomly oriented (parallely or perpendicularly) rela- 
tively to the external electric field, which is directed along x axis. In order 
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(a) (b) 

Fig. 2. a) An elementary square of the model with a vertical orientation, the dotted 
regions denote layers of the second phase; b) a lattice of the model, the small lines on 
the squares denote their orientations. 



for system to be isotropic the probabilities of the parallel and perpendicular 
orientations of squares must be equal or (what is the same) the concentrations 
of the squares with different orientations must be equal p\\= p^ — 1/2. 

This structure can appear, for example, on the small macroscopic scales, 
when a random medium is formed as a result of the stirring of the two-phase 
mixture. The corresponding averaged parallel and perpendicular (or serial) 
conductivities of squares a\ \ {x) and a±^ {x) are defined by the following formulas 

cr|| (x) — xai + (1 — x)g2 = cr+(l + 2ez), 

-^(-)=f- + ^) '=a+i-^. (11) 
\ai (J2 J 1 — Zez 

Thus we have obtained the hierarchical representation of random medium (in 
this case a two-level one). On the first level it consists from some regions (two 
different squares) of the random mixture of the two layered conducting phases 
with different conductivities cri and (T2 and arbitrary concentration. On the 
second level this medium is represented as a random parquet constructed from 
two such squares with different conductivities a\ | and a± , depending nontriv- 
ially on concentration of the initial conducting phases, and randomly distributed 
with the same probabilities pi — 1/2 (Fig. 2). This representation allows us to 
divide the averaging process into two steps (firstly averaging over each square 
and then averaging over the lattice of squares) and implement on the second 
step the exact formula (4). This can be considered as some modification of 
the FMSA approximation. As a result one obtains for the effective conductiv- 
ity of the introduced random two-phase model the following formula, which is 
applicable for arbitrary concentration 



1 + 2ez 



-,1/2 



1 - 2ez 



(12) 



This function has all necessary properties, satisfies equation (2) and coincides 
with the second form from IHUHI- 
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It is interesting to compare this formula with the known general formulas. 

(a) In case of small concentration of the first phase x <^ 1 one gets 

<^eff{x,{a})^a2(^l + ^^y (13) 

It follows from (13) that an addition of small part of the first higher conducting 
phase increases an effective conductivity of the system as it should be. 

(b) In the opposite case of small concentration of the second phase 1 — a; ^ 1 
one obtains 

ae/,(x,W)c.ai (^1-^^-1^), (14) 

i.e. an addition of the phase with smaller conductivity decreases (Jeff. It is 
worth to note that both these expressions for arbitrary values of the conduc- 
tivities CTi and (72 differ from equation (2) and coincide with it only in the 
weakly inhomogeneous case z <^ 1 . It must be not surprising because a form 
of the inclusions of the second phase in this model has completely different, 
layered, structure. In the low concentration expansion one can see again that 
the divergencies appear in the limit \z\ — s- 1. 

(c) In case of almost equal phase concentrations x — 1/2 + e, e^l one 
obtains 

aeff{e,{a}) ^ a+Vl - z^l + 2ez) . (15) 

The Keller - Dykhne formula (3) is reproduced for equal concentrations. 

One must note that at the same time the formula (12) does not satisfy the 
equation (6) except of the trivial case xi — X2. 

For a comparison of the different expressions of the effective conductivity the 
plots of /(e, z) in the EM approximation, in FMSA approximation and of the 
hierarchical model in FMSA-like approximation were constructed (see IHIE])- 
It follows from these plots that all derived above formulas for (Jeff , despite of 
their various functional forms, differ from each other weakly for z < 0, 8 due 
to very restrictive boundary conditions (6') and the exact Keller-Dykhne value. 
This range of z corresponds approximately to the ratio oijax ~ 10~^. For the 
smaller ratios a difference between these functions become distinguishable. 

Now let us consider in the more details the derived formulas for aeffi^, W}) 
in case when a2 — > 0{z — *■ 1). It is clear that for regularly inhomogeneous 
medium one can always construct such distribution of the conducting phase 
that cre//(e,l) will differ from zero for all 1/2 > e > —1/2. But in the case 
of randomly inhomogeneous medium the limit (72 — > is equivalent to the well 
known percolation problem [12,13]. In terms of z it corresponds to the limit 
z — > 1 and is also similar to the superconducting limit ai —^ oo . Strictly 
speaking, an implementation of the duality transformation (3) is not obvious in 
this case. However, if one supposes that the dual symmetry relation (3) fulfills 
in this limit too due to a continuity then it follows from (3) that 

t^e//(e)c^e//(-e) = 0. (16) 
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The relation (16) does not contradict to the known basic results of the percola- 
tion theory that de/ / (e) = for e < and (Jef /(e) for e > . Moreover 
one can show that in this case 



where a a is the odd part of cre//(e)- It means that a behaviour of aeff{e) in 
the percolation theory is completely determined by its odd part. It is known 
from the experimental and numerical results that in the percolation limit the 
effective conductivity a^f / have a nonanalytical behaviour near the percolation 
edge Xc — 1/2 (or at small e > ) 



where a critical exponent of the conductivity t is slightly above 1 and can be 
represented in the form t = I + S. Since the values of this exponent found by 
the numerical calculations are confined to be in the interval (1,10 - 1,4) [13], 
then S have to be small and belongs to the interval (0.1 - 0.4). It follows from 
(17) that the same behaviour must have cTq . It means that at small e there is 
some crossover on z under z —^ 1 from a regular (analytical) behaviour to a 
singular one. At the moment an exact form of this crossover is unknown. 

From the formulas obtained above, one gets always (Jeff "> in the limit 
a2 , except the region near a; = 1. It means that all these formulas, 
obtained in FMSA approximation, are not valid in the limit CT2 ^ 0. This is 
confirmed by the appearance of the divergencies in the limit z ^ 1. This fact 
is a consequence of the made approximation. For example, in case of the model 
of the layered squares this is due to the "closing" (or "locking") effect of the 
layered structure in the adopted approximation in the limit a2 0. In order 
to obtain a finite conductivity in this model above threshold concentration Xc 
one needs to take into account the correlations between adjacent squares. It 
is easy to show that near the threshold an effective conductivity is determined 
by random conducting clusters formed out of the crossing random layers from 
neighbouring elementary squares. As is well known, the mean size of these 
clusters diverges near the percolation threshold E| and for this reason the 
FMSA approximation cannot be applicable for the description of (Jeff in the 
region z ^ 1 and x < 1/2. It follows from our results that EM approximation 
overestimates aef / , whereas both other formulas underestimate it in this region. 
We hope to investigate the percolation limit in detail in the subsequent papers. 

Thus, though both formulae for the effective conductivity obtained above 
have the various functional forms, they satisfy all symmetries, including the 
dual symmetry and all necessary inequalities, and reproduce the general formu- 
lae for (Te// in the weakly inhomogeneous case. These results allow us to make 
a conjecture that aef / even of the two-phase randomly inhomogeneous systems 
may be a nonuniversal function and can depend on some details of the structure 
of the randomly inhomogeneous regions. The obtained formulae can be consid- 
ered as the regular ones, since they are applicable only for systems with ai ^ 0, 




(17) 



f^effie) ~ cri{x - XcY ^ CTie*, 



(18) 



7 



when there are no singularities connected with a percolation problem. The in- 
troduced composite method of the construction of the model random medium 
can be generalized on the heterophase systems with arbitrary number of phases 
and on the other ways of determination of the effective intermediate conducting 
boxes. It can be done for various types of boxes as well as for different numbers 
of the possible types of the boxes. It is clear that then one will have to use in- 
stead of (4) another formula. The constructed models and obtained expressions 
can be used for the modeUing and description of some real composite systems. 

The author thanks referees for useful remarks. This work was supported by 
RFBR grants # 2044.2003.2 and 02-02-16403. 
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